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In this paper, we investigate the solutions for a generalized fractional diffusion equation that 
extends some known diffusion equations by taking a spatial time-dependent diffusion coefficient and 
an external force into account, which subjects to the natural boundaries and the generic initial 
, condition. We obtain explicit analytical expressions for the probability distribution and study the 

relation between our solutions and those obtained within the maximum entropy principle by using 
Ch ' the Tsallis entropy. 
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I. INTRODUCTION 

Recently, anomalous diffusion equations have been extensively investigated due to the broadness of their 
physical applications. In fact, fractional diffusion equations and the non-linear fractional diffusion equations 
have been successfully applied to several physical situations such as percolation of gases through porous 
\jQ ' media [1], thin saturated regions in porous media [2], in the transport of fluid in porous media and in viscous 
' fingering [3], thin liquid films spreading under gravity [4], modeling of non-Markvian dynamical processes in 
protein folding [5], relaxation to equilibrium in system (such as polymer chains and membranes) with long 
temporal memory [6] , anomalous transport in disordered systems [7] , diffusion on fractals [8] , and the multi- 
1 physical transport in porous media, such as electroosmosis[9-I0]. Note that the physical systems mentioned 
above essentially concern anomalous diffusion of the correlated type (both sub and super-diffusion; see [11] 
. — ' an d references therein) or of the Levy type (see [12] and references therein). The anomalous correlated 
diffusion usually has a finite second moment (x 2 ) cx t a (a > 1, a = 1 and < a < 1 correspond to super- 
diffusion, normal diffusion and sub-diffusion, respectively; a = corresponds basically to localization). 
Due to the broadness of the problems involving anomalous diffusion, one needs to apply different kinds 
of theoretical approaches such as nonlinear Fokker-Planck equation (or modified porous media equation), 
fractional Fokker-Planck equation, Fokker-Planck equation with spatial dependent diffusion coefficient , and 
generalized Langevin equations. The properties concerning these equations have been intensively investigated 
[13-17] and the lattice Boltzmann method was used to get the numerical solutions for this equations which 
govern the multi-physical transfort in porous media. 

In order to cover the above situations, we employ a spatial time-dependent diffusion coefficient, in other 
words, our work is aimed at the investigation of solutions for a fractional diffusion equation taking a spatial 
time-dependence on the diffusion coefficient and an external force (drift) into account. More precisely,we 
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focus our attention on the following equation: 

J^, *) = f o dt '^ D ^ 1 - <£^r *)H - ^m*M*, 0}, (i) 

where < 7 < l,p, v € R, the diffusion coefficient is given by D(x,t) = D(i)\x\~ e , which is a spatial time- 
dependent diffusion coefficient, and F(x) = — dV J^ is an external force (drift) associated with the potential 
V(x). Here, we use the Caputo operator [18] for the fractional derivative with respect to time t, and we 
work with the positive spatial variable x. Later on, we will extend the results to the entire real x-axis by 
the use of symmetry (in other words, we are working with d/d\x\ and d^ 1-1 /dlx^ -1 ). Also, we employ, 
in general, the initial condition p(x,0) = p{x) (p(x) is a given function), and the boundary condition 
p{x — >±oo,t) — > 0. For Eq.(l), one can prove that f_°° dxp(x,t) is time independent (hence, if p(x,t) is 
normalized at t = 0,it will remain so forever). Indeed, if we write Eq.(l) as = —d x J, and, for simplicity, 
assume the boundary conditions J(±oo) = 0, it can be shown that dxp(x, t) is a constant of motion (see 
[19] and references therein). Note that when (p,j, v) = (2, 1, 1), Eq.(l) recovers the standard Fokker-Planck 
equation in the presence of a drift taking memory effects into account. The particular case F(x) = (no drift) 
and D(x,t) = DS(t) with (p,j) = (2,1) has been considered by spohn [20]. The case D(x,t) = DS(t)\x\~ e 
with (p, v) = (2, 1) and the case D(x, t) — D{t) with (p, v) = (2, 1) have been investigated in [21] and [22], 
respectively. 

Explicit solutions play an important role in analyzing physical situations, since they contain, in principle, 
precise information about the system. In particular, they can be used as an useful guide to control the 
accuracy of numerical solutions. For these reasons, we dedicated to this work to investigate the solutions 
for Eq.(l) in some particular situations. In all the particular cases, Eq.(l) satisfies the initial condition 
p(x,0) = p(x) (p(x) is a given function), and the boundary condition p(±oo,t) = 0. The remainder of this 
paper goes as follow. In Sec. 2, we obtain the exact solutions for the special cases. In Sec. 3, we present our 
conclusions. 



II. EXACT SOLUTIONS FOR DIFFERENT CASE 



In this section, we start our discussion by considering Eq.(l) in the absence of external force with 
D(t) = Dt a ^ 1 /T{a) ( D(t) = D6(t)),(p,v) = (2, 1) and 7, 6 arbitrary. For this case, Eq.(l) reads 

^p(*,*) = ^^(t-0^{N-^p(*,*)}. (2) 

Here, we use the Caputo operator [18] for the fractional derivative with respect to time t. By employing the 
Laplace transform in Eq.(2), we obtain 

D(s)^-{\x\- e ^p(x, s)} - s^Pix, s) = s^pix, 0), (3) 

where p(x, s) = J^{p(x, t)}, D(s) = JC{D(t)}, and jSf {/(£)} = J °° dte~ st f(t) denotes the Laplace transform 
of the function /. This equation can be solved by Green function method [23]. By substituting 

p(x,s)= / dx'Q{x — x', s)p(x'), (4) 
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into Eq.(3) which yields 

^)^{\x\- e ^g{x lS )}-s^g{x lS ) = - s ^- 1 g(x,o). (5) 

where Q(x,t) subjects to the initial condition Q(x, 0) = S(x) and the boundary condition Q(±oo,t) = 0. 

In order to solve Eq.(5), it is convenient to perform the transform [24] 

y = A( S )x\ g(x lS )=y 5 Z{y) (6) 

to translate Eq.(5) into the second-order Bessel equation as 

2 d 2 Z dZ 2 2 y 2 - A y i 

y W + V &H (A + V )Z{V) = -— S{i W) } " } (?) 

with parameter A 2 under the following conditions: 

2 + 6> N l r sT ,i , l + 6>. l + 6> 

Since Eq.(5) should fit the boundary condition Q(±oo,t) = 0, i.e. C/(±oo,s) = 0, we get the solution of 
Eq.(5) 

g(x,s) = C(s)y 5 K x (y), (9) 

where K n (x) is the modified Bessel function of second kind; and C(s) can be determined by the normalization 
of g(x, t), i.e. Jp 00 dxg(x,s) = j^. After some calculations, we obtain 

r (2Te)- s 2 + 61 W ^ v 2 + 6» v £)(s) ; 1 1 ; K ' 

where, we used the formula 

' dy • y^{ay) = 2-*«— * r(i±f^ )r (i±^). (11) 
Case 1. D(t) = D8(t), i.e. D{s) = D. 

Since K\{x) = 2[tI(-a/2,i)(a/2,i)]. we can get the Laplace inverse of g(x,s) by applying the 

property of the Laplace inverse of Fox function, which yields 

a(rt)- 2 + 9 ( 1 \^h 2 or M 2+e la-sfr.T) i (m 
y[x > l) ~ 2 r(^) [ (2 + 6) 2 Dr» ) 1 2l (2 + fl) 2 5n l WD,(^,D 1 ' 1 j 

where if™ g ™M(b 1 I 's^''"'(b g * ] is the FOX function [25]. Thus, we can find the solution by substituting 
Eq.(10) into Eq.(4), which yields 

p(xt) - 2 + 9 ( 1 ^ r + °°dx'p(x')H 2 or k~^| 2+e .d-.fr, t) , (13) 
2 r(2^) t (2 + 0) 2 J D^ j 7.^ dx P^)*i 2l( 2 + fl)2 i?4 7l(o,i),(^.i) J - ^ 

In fig.l, we show the behavior of the above equation by considering typical values of 7 and 9 with /5(x) = (5(a?). 
At this point, it is interesting to analyze the asymptotic behavior of Eq.(13). For simplicity, we consider 
p{x) = S(x), so p{x,t) = g(x,t)] and the asymptotic behavior of p(x,t) is 

2 + 9 1 7 _i 1 1 -,-1 

0(X t ) ~ (2 — ^)~2^(2 + 0)(2- 7 ) 2 C ) (2 + 0)(2- 7 ) | X |5=7 

^ ' J 2r(l/(2 + 6')) 1 7j 7 1 (2 + 0)2£>*7 j 11 

\x\ 2+e _L_ 

x eap(-(2 - 7 )7 2 -' ( (2 I e y mi ) 2 ' T )■ ( 14 ) 
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and 



In this direction, Eq.(14) can be considered as an extension of the asymptotic behavior of homogeneous and 

isotropic random walk models [26] . 

Case 2. D(t) = , i.e.D(s) = D S - a . 

By using the same method as in case 1, we obtain 

C ixt )- 2 + 9 ( 1 )^H 2 or M 2+<) ,(1-^.7+^ , , 

y[x ' t} ~ 2 r(^) l (2 + 0) 2 ^+« j 1 2[ (2 + 0)2 J Di7+« l (O4),(i±|,i) J ' [W) 

, \ 2 + , 1 ,i f + °° , , NTT , nr |x-a:'| 2 + 8 | (i-i+e, 7+a ) 1 

^■ t) = ^(i5) ( (2 + fl )^ ) ^/, BO ^H^J wS') 1 ' (16) 

Let us go back to Eq.(l), and consider the external force F(x) oc x\x\ a ~ l , D(x,t) = DS(t)\x\~ e and 
/i = 2, v = 1. In this case , analytical solution can not easily be obtained for a generic a, 0. However, for 
6 7^ 0, and a + 6* + 1 = 0. By following the same procedure as in the above case, an exact solution can be 
obtained and it is given by 

o(xt) - 2 + 6 ( 1 )^ g 2 Or N 2 ^ ,(1-^,7) , (17) 

where, for simplicity, we are considering the initial condition p(x,0) = 5(x), and the external force (drift) 
F(x) = K.x a . The second moment is given by < x 2 >oc t^+e, which scales with the exponent and clearly 
depends only on 7 and 6. So, when < 1, = 1 and > 1, the system is sub-diffusion, normal diffusion and 
supper-diffusion respectively. 

The presence of the external force in Eq.(l) is now changed into Fix) — —k\x + fex -1- e . In order 
to obtain the solution of Eq.(l), we expand p(x,t) in terms of the eigenf unctions, i.e. we employ p(x, t) = 
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St^Lo ( t >n 0*0 with ip n (x) determined by the spatial equation 



- \nMx) = D ^ x ~ e ^ n ^ + §i {{klX + k ^~ 1 ~ e )M^)} (is) 

and <fi n (t) determined by the time equation 

— Mt) = ~K4>n(t)- (19) 
The solution for the time equation is given by in terms of the Mittag-Lcffter function 

<f>(t) oc £ 7 (-A„^). (20) 
In order to get the solution for Eq.(18), we perform the transform 

iMx) = e-VZ(y), y = Ax v (21) 
to translate Eq.(18) into the associated Laguerre equation 

yZ"{y) + (a + l- y)Z'{y) + nZ(y) = (22) 
with parameter a under the following condition: 

fcg 1 f\ I T 

a = — — , v = 2 + 9,A= 2 -—,5=- 1— ,A„ = (2 + 6)nki. (23) 

2 + 6* (2 + 9)D' (2 + 9)D' n K J y ' 

Then, using the Green function methods and after some calculations, it is possible to show that 

/>oo 

p(x,t) = / dx p(x )G(x,x ,t), (24) 



ki fc 2+p , k 2 __^_| r |2+f< ^2, (2 + 0)T(n + 1) , , 

5(^X0,0 = (— 1— )<™|x|^e E r , ! l+D E,(-^n 

y^^Ujiy n=Q L[ - n+ (2+8)D> 

where Ln\x) is the associated Laguerre polynomial and it is the solution for Eq.(22). Here, we used the 
formula 

' dyy<*e-yLl«\y)Lt\y) = r( " + ° + 1) . (26) 

o 1+1 

Notice that Eq.(25) contains the usual Ornstein-Uhlenbeck process [27] and the usual Raylcigh process [28] 
as particular cases and it extends the results obtained in [21]. In this context, in the presence of an constant 
absorbent force, i.e. ap(x, t), we can obtain the solution for Eq.(l) only need to change the argument of the 
Mittag-Leffler function present in Eq.(20) to A n + &■ 

Let us now discuss Eq.(l) by considering a mixing between the spatial and time fractional derivatives. 
For simplicity, we consider Eq.(l) in the absence of the external force with (y, 6) = (1, 0) and 7, p arbitraries. 



6 



Applying the Fourier and Laplace transform to Eq.(l) and employing the Riez representation for the spatial 
fractional derivatives, we have 

s 7 ,5(/c, s) - s r "V(fc, 0) = -Dtflkfpfc s), (27) 

where p(k,t) — J-{p{x, t)} = p(x, t)e~ %kx dx , so p(k,0) = 1. Here, we consider the diffusion coefficient 
given by D(t) = t^o) » *' e ' ^( s ) = Ds~ a and 0) = By using the inverse of Laplace transform, we 

obtain 

p(M)=£ 7+ a,l(-#|fcR 7+Q ) 



= H\ ^Dim^Xlio^- (28) 

This solution recovers the usual one for (/x, 7) = (2, 1) and for ix 7^ 2 it extends the results found in [29]. In 
order to perform the inverse of Fourier transform, we employ the procedure presented in [30] . Then we can 
obtain 

n f„ +\ _ \ tt2 lr N i(l-l/M,l/M),(l-(7+")/M,(7+a)/M)l /%\ 

^ ' j " 2|V5F(i^)- 3[ 2( ^ 7+q) J ( 0,1/2),(1-1/,,1/,),(1A1/2) J' ^ 

The stationary solution that emerges from this process is a Levy distribution. 

Now, we consider a particular case of Eq.(l) for 7 = 1 and nonzero values of p and 6, and consider a 
linear drift, i.e. F(x) = —K.X. For simplicity, we employ D(t) = D6(t) and the initial condition p{x, 0) = S(x), 
then Eq.(l) yields to 

8 8 8^~ 1 8 

di p{ - x ^ =D d- x { ^ ^w^)n+^{^p(M)}- (30) 

Let us investigated time dependent solutions for Eq.(30). We use similarity methods to reduce Eq.(30) 
to ordinary differential equations. The explicit form for these ordinary differential equations depends on 
the boundary conditions or restrictions in the form of conservation laws. In this direction, we restrict our 
analysis to find solution that can be expressed as a scaled function of the type [31] 

p(x ,t) = ^p(z), z=^ fy (31) 

Inserting Eq.(31) into Eq.(30), we obtain 

. (bit) K s d . D 8 . g d^ 1 , U1 

- W + m ] 8- z [zp{z)] = m^^ [z d^ p(zn (32) 



By choosing the ansatz 

<b{t) K, kD 



(33) 



cb{t) 2 cb(t) <b{t) 9 +v+"' 

where k is an arbitrary constant which can be determined by the normalization condition. By solving 
Eq.(33), we have that 

(bit) = [U( Q )f +^+,-1^(6+^-1)^ + Dk _ ^(fl+^-Dtt)]^^ (34) 

K, 
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By substituting Eq.(33) into Eq.(32), we obtain 

Ft cl^-i ft 

a- z [z a^^ = - k ai [zm] - (35) 

Then, we perform an integration and the result is 

z- 9 ^~p{z) v = -kzp(z)+C, (36) 
where C is another arbitrary constant. Also, we use the following generic result [32]: 

D s x [x a (a + bxf] = o* r[ "+ 1] s°-*(a + bxf~ 5 (37) 
1 [a + 1 — o\ 

with D & x = d s / dx & and 8 = a + /3 + l. By defining g(x) = x^ (a + bx)" and A = a(l — ^) — 8, and rearranging 
the indices, Eq.(37) can be rewritten as follows: 

d 5 Mx)Y = aa r[a + i-V Ag(a:) ' (38) 

For this case, we consider the ansatz p(z) = Afz^(l + bz)». By using the property of Eq.(38) in Eq.(36) 
and ,for simplicity, choosing C — 0, we find 

(2- M )(/i + 0) 
a 1-2/1-0 ' 

2-n 

v = 



1 + fX + i 



In this case, we have 

= ^) [ (i + 6z)(^)aw) ]l """ e ' (40) 

p/ Q\ f.t-\-0-\-l 

where <p(t) is given above, Af = [~fc r(a+i) ] 1 ~ 2m ~ 8 an d & is an arbitrary constant (to be taken, later on, as 
±1 according to the specific solutions that are studied). Several regions can be analyzed. For simplicity, we 
illustrate two of them: — oo < fj, < — 1 — 8 with 9 > 0, and < fj, < 1/2 with < 8 < 1/2 — \x. Let us start by 
considering the region — oo < \i < — 1 — 6. Without loss of generality, we choose b = — 1. The normalization 
condition implies (see Fig. 2) 



So 



N \l l-Tn uTT TS^ ] r= ^dz= 1. (41) 

7_1 (1 + 5 Z )(1-M)(1+M+0) J V > 



T\l-u-9] , , 

AT = t — J (491 

J r, r r ^ 2 +A'g-2e-2^ i r r l-^+^ 2 +e 2 +2^e i • V ^ 
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FIG. 2: Behavior of <j>(t)p(x,t) versus x/<f>(t), which illustrates Eq.(40) with typical values for p and 9 satisfying 
()</!<— I — 9 and 9 > 0. We notice that the distribution vanishes at the abcissa equal ±1, and remains zero outside 
of this interval. 
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FIG. 3: Behavior of <f)(t)p(x,t) versus x/(f>(t), which illustrates Eq.(40) with typical values for /j, and 9 satisfying 
< /J. < 1/2 and < 9 < 1/2 — p. We notice that the distribution vanishes at the infinity. 



Let us now analyze the region < fx < 1/2 with < 6 < 1/2 — pi. Again without the loss of generality, we 
choose b = 1. The normalization condition implies (sec Fig. 3) 



r 



i+e-M 



Let us finally mention a connection between the results obtained here and the solutions that arise from 
the optimization of the nonextensive entropy [33] . These distributions do not coincide for arbitrary value of 
x. However, the comparison of the \x\ ~ > oo asymptotic behaviors enables us to identify the type of Tails. 
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By identifying the behavior exhibited in Eq.(40) with the asymptotic behaviors l/|x| 2 '^ 9 that appears in 
[33] for the entropic problem, we obtain 

q = TTJ+8' (44) 

This relation recovers the situation for 6 = . 



III. SUMMARY AND CONCLUSIONS 



We have analyzed the generalized fractional diffusion equations by considering an external force F(x) oc 
x|x| a_1 and a spatial time-dependent diffusion coefficient D(x,t) = D{t)\x\~ e . By using Laplace transform, 
Fourier transform, the Green function method and normalized scaled function we can find the explicit 
solutions p(x,t) which subjects to the natural boundary condition p(±oo,i) = and the initial condition 
p(x,t) = p(x). In a word, we have extended the results previously obtained by the other authors by taking 
an external force and a spatial time-dependent diffusion coefficient into account. Wc have also discussed 
the connection with nonextensive statistics, providing the relation between our solutions and those obtained 
within the maximum entropy principle by using the Tsallis entropy. Finally, we expect that the results 
obtained here may be useful to the discission of the anomalous diffusion systems where fractional diffusion 
equations play an important role. 

Acknowledgments: The authors are grateful to the anonymous referees for useful comments and 
suggestions. 
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